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Abstract. Resonant enhancement of high harmonic generation can be obtained 
in plasmas containing ions with strong radiative transitions resonant with 

harmonic orders. The mechanism for this enhancement is still debated. We 
perform the first temporal characterization of the attosecond emission from a tin 
plasma under near-resonant conditions for two different resonance detunings. We 
show that the resonance considerably changes the relative phase of neighbouring 
harmonics. For very small detunings, their phase locking may even be lost, 
evidencing strong phase distortions in the emission process and a modified 
attosecond structure. These features arc well reproduced by our simulations, 
allowing their interpretation in terms of the phase of the recombination dipole 
moment. 
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1. Introduction 

The strong-field process of high harmonic generation (hhg) provides a source of 
coherent femtosecond and attosecond pulses in the extreme ultraviolet (xuv) with 
a wide variety of applications [T. It has allowed the study of photoionization of 
atoms and molecules with attosecond resolution [21 [3], revealing unexpected features, 
in particular near resonances [H |5] . One would expect that such resonances also affect 
the HHG process itself, since the final step of the conventional semi-classical model [6l[7j 
is photo-recombination, i.e. the inverse process of photo-ionization [511^110]. However, 
few experiments have shown evidence of resonance effects [TTl [TH [131 [HI [H] i until 
recently when hhg was demonstrated from ablation plasmas produced on solid targets 
|16j . This technique allows using a wide variety of ionic targets for hhg, including some 
with very strong radiative transitions. Enhancement of single harmonics resonant with 
such transitions has resulted in conversion efficiencies of ~ 10"'* [12l[l8], far exceeding 
the typical value of ~ 10^^ obtained in neutral gases TO]. 

This resonance-enhancement is currently a hot topic of both applied and 
fundamental interest. However, so far, only the spectral intensity of high harmonics 
from resonant plasmas has been experimentally characterized and the important 
question whether they offer a route to intense XUV femtosecond or attosecond 
pulses remains unanswered. More fundamentally, theories describing the still scarcely 
explored role of resonances, and thus atomic structure, in strong-field processes 
[ini[ni[21l23IM[151[lSllS], would greatly benefit from advanced experimental 
characterization of the emission as a benchmark. As yet, it was shown that the 
resonance-enhancement effect vanishes when shifting the harmonic energy off the 
resonance [T71 [T51 [57] • Furthermore, the xuv yield was found to decrease with 
increasing ellipticity of the driving laser |17j. suggesting electron recoUision as a crucial 
part of the enhancement mechanism. The ellipticity dependence is, however, much 
slower than for rare gas media [28 . 

In this article, we report measurements of both the femtosecond and attosecond 
structure of hhg from tin ablation plasma, where harmonic 17 (HIT) of an « 800-nm 
laser is tuned into resonance with the 4d**'5s^5p ^P3/2 4d^5s^5p'^ {^D)^D^^2 
transition in Sn''" with energy Cros = 26.27 eV. We find a femtosecond envelope 
that does not correspond to a slowly-decaying plasma emission but is indeed slightly 
shorter than the driving laser pulse. However, the harmonic relative phase, which 
governs the attosecond structure, is significantly perturbed by the resonance. In 
resonant conditions, the phase locking between the resonant and the neighbouring 
orders is lost, i.e. their relative phase varies significantly within the harmonic spectral 
width. We develop a model based on the numerical solution of the time-dependent 
Schrodinger equation (TDSE) that reproduces the experimentally observed strong 
phase perturbations, both on the attosecond and the femtosecond scales. Furthermore, 
we clearly link these effects to the resonance by means of an analytically calculated 
recombination dipole matrix element. 

2. Experimental study 

Our experimental setup is similar to the one in I29j . The 20-Hz, 50-mJ Ti:sapphire 
laser is split into two beams before compression: a w 300-ps pre-pulsc is loosely 
focused to ~ 10*° W cm~^ intensity onto the solid tin target at normal incidence, 
creating an ablation plume of < 1 mm length. The second beam is delayed by 
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Figure 1. Spectra obtained by photoionizing Ar gas with high harmonics 
generated in tin plasma for the 'detuned' (a) and 'resonant' (b) cases, in presence 
of the IR probe beam. We corrected for the ionization cross-section such that 
the quantity shown is proportional to the xuv spectral intensity, except for the 
presence of SBs, that only appear in the photoelectron spectrum. The 0.16 eV 
wide red line marks the resonance and the dashed line marks the exact position 
of H17. The inset in panel (a) shows the theoretical photoionization cross section 
of Sn+ 27 , on the same energy scale as the HHG spectrum. 

80 ns and then compressed to tjr = 55 fs duration, before being focused (//50) to 
/J^™^ < 1 X lO^^'W crn"^ (estimated by measuring the attochirp of the harmonic 
emission from an argon gas jet in the same experimental conditions 30J) into the 
plasma plume, propagating at a distance of 300 jim parallel to the target surface. 
The annular far-field profile of the HHG driving beam allows blocking it with an iris 
« 1 m downstream, while the harmonics pass on-axis. The generated harmonics are 
detected by photo-ionizing argon atoms in a magnetic bottle electron spectrometer 
(mbes). Here, a weak ir probe beam (/{^°''° ~ 10^^ Wcm~^) can be overlaid, leading 
to two-photon xuvzbiR ionization. Thus, between the odd-order harmonic peaks, 
even-order sidebands (SBs) appear in the photo-electron spectrum |31| . which serve 
as the nonlinear signal for the temporal characterization. In these very challenging 
experiments, the signal-to-noise ratio was limited by the relatively small number of 
pre-pulses (few thousands) that the tin target could withstand before deterioration 
and drop in harmonic signal [29) . 

Two separate measurement runs with the laser central wavelengths Al — 793 nm 
and 796 nm were made, corresponding to detunings of H17 from the resonant 
wavelength Arcs — c/i/crcs of jAics — Air/17|/Ai.os = 7.5 x 10^"^ and 1.13 x 10~^, 
respectively. While this difference is small, it shifts the H17 photon energy by 0.1 eV, 
which is significant compared to the 0.2 -eV typical width of a harmonic line and the 
« 0.16 -eV resonance width [2^. The relative positions of H17 and the resonance for 
the two driving wavelengths are marked in figure [T] for Al = 796 nm, H17 is closer 
to the exact resonance than for 793 nm, and for the sake of brevity, we will in the 
following use the oversimplified terms of 'resonant' and 'detuned' case, respectively. 
Note that we expect the resonance to play a role in the emission of H17 in both 
cases. We do not consider the AC Stark shift of eros which is expected to be small at 
our laser intensities [57^. A possible small blue-shift of Al due to a time-dependent 
electron density during the driving pulse is also neglected. No such blue-shift, or 
other significant phase distortions have been observed in spider measurements of the 
driving laser pulses after propagation through the plasma plume. Note finally the array 
of 4d O 5p transitions around 26 eV (see inset in figure [T^) . The above-mentioned 
transition to the (^D)^D5/2 autoionizing (Al) state is expected to play the dominant 
role because of its gf -value exceeding those of the others by a factor > 3.5 p7] . 
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Figure 2. Temporal characterization of the harmonic emission from tin plasma 
plume. (a)/(b) Measured data (circles) with 5-point Savitzky-Golay smoothing, 
from low resolution XUV-IR cross-correlations for SBs 14/18 in the 'resonant' case. 
The gray shaded area is an intensity auto-correlation of the IR laser pulse, (c) 
Attosecond pulse profile (gray shaded area), in the 'detuned' case, reconstructed 
from the RABBIT scan shown in figure |3j3,c, taking into account harmonics 13 to 
19. The dashed line is the IR driving field modulus. 

The distinct difference in resonance conditions is evident in tfie hhg spectra shown 
in figure [l] The laser focusing geometry, pulse energies and durations were kept the 
same in both runs, which ensures equal intensities for the pre-pulse and driving pulse, 
as well as equal phase matching conditions, assuming similar spatial beam quality. 
However, H17 shows significantly stronger enhancement in the 'resonant' case, while 
H13 and H15 are of essentially the same intensities in both cases. In the 'detuned' 
case, varying the generation conditions and in particular increasing the hhg driving 
intensity did not allow us to increase the H17 relative intensity. We thus attribute the 
difference to the driving laser wavelength. 

The enhancement observed in figure [l|^b) is smaller than that reported in [17] . 
This may be due to the different macroscopic plasma medium conditions, caused by 
our relatively large (~ 300 jim) distance of the driving beam from the target (so as to 
avoid high electron density gradients that diffract the IR light, hampering the blocking 
of the driving beam) , and smaller pre-pulse spot size resulting in a shorter generating 
medium. 

2.1. Femtosecond characterization 

In a first step, low-resolution (3.3 fs delay steps) xuv-iR cross-correlation scans 
were performed to measure the xuv pulse envelopes. The SB intensity /sb(''') = 
/_oo ^x\Jv{t)f[hR{t — t)] dt provides a cross-correlation of the xuv and the IR pulses, 
with temporal intensities Ix\jv(t) £^nd IiB.(t), respectively [31| . The function /[/jr] 
is well approximated by a power law cx /ir" at low IR intensity; in our experimental 
conditions, we find a ~ 0.8 [35] • Figures [2j a) and (b) show the intensities of the 
SBs 14 and 18 as a function of the Xuv-IR delay t in the 'resonant' case. Given 
the high intensity contrast of H13 vs. HIS and H17 vs. H19, these two SBs 
essentially compare the durations of H13 and H17, i.e. a non-resonant with a resonant 
harmonic. Within the experimental precision, both cross-correlation traces have the 
same shape and width (tsb — (80 ± 10) fs FWHM) as the intensity auto-correlation 
of the 55-fs laser pulse (measured using a standard second-order auto-correlator). 
Data from the 'detuned' case show a very similar behaviour ^32. . The harmonic 
duration, txuVi can be estimated by assuming Gaussian intensity profiles, such that 
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Txuv = (tsb^ - TiRVa)i/2 [31]. We then find txuv = (50 ± 15) fs for both H13 
and H17. For a 55-fs driving pulse, this is fairly long, especially for harmonics close 
to the cutoff. A likely reason is that hhg was saturated before the driving pulse 
peak. Indeed, the observed hhg cutoff at H19 suggest an effective driving intensity 
of > 6 X 10^^ W cm"^, which is already higher than the barrier suppression intensity 
[33] for the ionization of Sn+ (/p = 14.6 eV). 

2.2. Attosecond characterization 

In order to access the attosecond sub-structure of the hhg emission, the xuv-iR cross- 
correlation scans have been repeated varying r with sub-laser-cycle resolution (0.15-fs 
delay steps) over sa 7 f s around the peak overlap. The IR probe beam intensity has 
been reduced by a factor 2 to completely suppress higher than first-order SBs [34) . 
This is an implementation of the rabbit technique [33]: the high resolution in r 
allows observing oscillations of the SB intensities with twice the IR laser frequency, 
wl, as a result of interference between the two quantum paths involving absorption 
of the neighbouring harmonic photons and either absorption or emission of an IR 
photon. From the phase of this oscillation in each SB g, the relative phase, A(/5q, of 
two neighbouring harmonic orders [q — 1) and {q + 1) can be extracted. This yields 
an approximation to the xuv group delay (GD) dLp/duj\^^^ A(/3^/2wl, where wl is 
the laser frequency. 

More precisely, since we are not dealing with discrete harmonic lines, one can write 
the r-dependent sideband intensity spectrum as SB(uj,t) = \A{uj + wl) exp[ia;LT] -|- 
A(a; — wl) exp[— I^ltJ^, where A{u!) is the harmonic complex spectral amplitude. The 
integration over uj of the SB line, as usually done in experiments (the convolution with 
the probe-iR spectrum (O.OSeV FWHM) and with the spectrometer resolution (0.25eV 
FWHM) in any case washes out the w-dependence of the 2cjL-oscillation within the 
SB line), thus gives access to a mean value of the harmonic phase difference if the 
latter does not vary significantly over the harmonic spectral width: the two harmonic 
orders are then said to be locked in phase. A significant variation would result in a 
disappearance of the 2ajL-oscillation: this has never been observed in any media, be 
it atomic gases [SS] 130] , molecular gases [Ml I3Z] or non-resonant plasmas [59] . Such a 
disappearance would correspond to significant variations of the individual attosecond 
pulses throughout the emitted pulse train [38l [39] . 

Figures [sja) and (b) show the squared Fourier transform of the r-dependent 
integrated SB intensities from rabbit scans for the 'detuned' and 'resonant' case, 
respectively. In the 'detuned' case, the 2wL-oscillation is detected with reasonable 
contrast on all observable SBs. In the 'resonant' case, we find a striking peculiarity 
of the two SBs (16 and 18) surrounding the resonant HIT: though present in the 
photoelectron spectrum and of similar amplitude as in the 'detuned' case (see figure 
[T]), they do not oscillate at 2ajL- This oscillation is detected with good contrast on 
all other SBs to which the resonant HIT does not contribute, demonstrating that 
interferometric stability was sufficient. The strong intensity difference between H17 
and HIS, that leads to a reduced interference contrast, is the same between H15 
and II13 and does not prevent observation of the 2a;L-oscillations of SB14. We thus 
attribute this suppression to the closer proximity of HI 7 to the exact resonance. 

From the SBs with detectable 2wL-oscillation, the GDs shown in figure |3jc) are 
obtained. For SB12 and SB14, the measured GDs in the 'resonant' and 'detuned' 
cases agree well within the experimental precision, underlining that the experimental 
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Figure 3. Results of rabbit measurements. Squared modulus of the Fourier- 
transformed oscillating SB signals in the 'detuned' (a) and 'resonant' (b) cases. 
The phase at the peaks near 2ajL gives the GDs, shown in (c) for the 'detuned' (♦) 
and 'resonant' cases (□). Errors are the standard deviation of the phase within 
the FWHM of the 2ajL-peak. Red dots and dashed line show GDs extracted from 
numerical simulations for the 'detuned' case (see text) under a driving intensity 
of 6.5xl0i3 W cm . Solid line shows an SFA-calculation |30| for this intensity 
and the ionization potential of Sn+. 

conditions have been equal, except for the small laser wavelength shift that tunes 
H17 into resonance. This also shows that, despite the reduced signal-to-noise ratio, 
leading sometimes to small shifts of the 2wL-peak positions, the measured phases are 
reliable. The GDs at SB12 and SB14 also agree well with those predicted by the 
strong field approximation (SFA) for the short trajectory class and a driving intensity 
of 6.5 X 10^"^ Wcm^^ (solid line in figure [sjl^c)). For SB16 and SB18, we can only extract 
GDs for the 'detuned' case. For SB16, we find a striking increase of the GD by 350 as 
as compared to the SFA prediction. As we will show, this strong phase distortion is a 
result of the nearby resonance. For SB18, the measured GD is much closer to the SFA 
prediction again, even though a big instability for this specific SB was observed in the 
different data sets. Finally, for SB20 we only obtain a good SB oscillation contrast in 
the 'resonant' case and for a single rabbit scan. The extracted GD carries a rather 
large error due to the low signal-to-noise ratio and may not be reliable. We will thus 
not discuss it further and show it only for completeness. 

In figure [2|c) , we show the average attosecond pulses in the train reconstructed 
through a Fourier transform from the measured GDs and intensities of harmonics 13 
to 19 in the 'detuned' case. The rather big attochirp together with the phase jump 
between H15 and H17 (corresponding to the GD shift of SB16) distorts the pulse 
shape to a double peak per half cycle. 

3. Theoretical study and discussion 

In order to interpret our experimental observations, we made simulations based on 
the numerical solution of the 3D-TDSE for a single active electron in a model ionic 
potential, modulated by a strong laser field. The potential (similar to the one 
shown as a black line in figure |4j see details in [201 HZ]) has a quasi-stationary 
state modeling the (^D)^D5/2 AI state that we expect to dominate the resonance- 
enhancement effect. Namely, we reproduce the ionization and the AI state energy, the 
width and (less accurately) the oscillator strength of the ground- AI state transition 
in real Sn+. The laser pulse has a cos^-intensity-envelope with a duration of 55 fs 
FWHM. The computed time dependent expectation value of the dipole acceleration 
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Figure 4. Illustration of the role of the AI state in resonant HHG. See text for 
details. 

gives the harmonic spectrum in amphtude and phase. 

In order to partly simulate the macroscopic effects, and in particular to select 
the shortest quantum path contribution, we perform a coherent summation of the 
numerical spectra over the laser intensity (steps of 2x10^^ W/cm^), thus calculating 
an on-axis XUV signal from a thin target [40l|41]. In this signal, the shortest trajectory 
contribution dominates |40) because phases of longer trajectories contributions vary 
rapidly with the laser intensity [H] . 

As it was suggested in [5D] and further studied in [571 HI] ; the role of the resonance 
with the AI state can be understood within the four-step HHG model illustrated in 
figure [4] After (tunnel) ionization (step 1) and acceleration in the laser field (step 2), 
the electron can recombine to the ground state and emit xuv (step 3* in the three- 
step model [7]); however, if the energy of the returning electron is close to that 
of the AI state, the electron can be captured such that the system finds itself in the 
AI state (step 3) and only then relaxes to the ground state, emitting XUV (step 4). 
Thus, the phase of the resonant harmonic is defined by the phase accumulated during 
the free-electronic motion and during recapturing of the electron. The calculation 
results presented below show that the latter is very sensitive to the detuning from the 
resonance. 

Group delays were found from the numerical spectrum by taking the phases at 
the exact harmonic frequencies, {2q± l)a;L and calculating the phase differences Aipq 
and thus the GDs A</5g/2a;L- The robustness of the extraction of GDs from the TDSE 
solutions has been verified by comparing to two additional methods: (ii) After selecting 
in the spectral domain only two neighbouring harmonics, (2g± 1), we make an inverse 
Fourier-transform back to the time domain, where we obtain a train of attosecond 
pulses. The temporal position of the first attosecond pulse peak after the driving 
pulse maximum gives the GD at SB q. (iii) Selecting in the spectral domain only single 
harmonic peaks, inverse-Fourier-transforming to the time domain and multiplying by 
ex.pliqwi,t], we find the q-th harmonic complex slowly varying envelope, aq{t). The 
harmonic temporal intensity and phase are then |ag(<)p and arg[aq(t)], respectively. 
We then calculate the GDs from the time averaged relative phases of neighbouring 
harmonics. The results of these three methods agree to within ±60 as. 

In figure [3]:, we show the so obtained GDs in the 'detuned' case for a laser peak 
intensity of Iq — 6.5 x 10^^ W cm~^. Very satisfactory agreement is obtained with 
the measured values. In particular at SB16, the obtained GD shift of « 500 as with 
respect to the SFA reproduces well the experimental one. This shift is quite robust in 
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Figure 5. Phase and squared modulus of the recombination dipole matrix 
element (x)c— g as a function of the emitted XUV photon energy (i.e. the difference 
of the continuum- and ground-state energy). The 0.2 eV wide blue and yellow 
lines mark the spectral position and approximate width of H17 for the detuned 
and resonant case, respectively. H15 and H19 are around 23.4 eV and 29.6 eV, 
respectively. 



the simulations for different peak intensities, Iq, as well as in all our measurements. 
In contrast, the GD at SB18 turns out to be very sensitive to Iq in the simulations 
and for some values, large shifts are obtained. We have also observed such shifts for 
SB 18 in a few measurements (not shown). 

The contrast of the 2a;L-oscillations of the different SBs was evaluated by 
integrating the above expression of SB{uj, r) over the spectral width. In the 'resonant' 
case, the contrast values for SB12 to SB20 are respectively: 0.69, 0.62, 0.08, 0.13, 0.45. 
In qualitative agreement with the experiments, the calculated contrast for SBs 16 and 
18 is strongly reduced. Moreover, it is significantly lower in the resonant than in the 
detuned case, while being approximately the same for the other sidebands. 

The above results can be understood by deriving analytically the wavefunctions 
(both in the discrete and continuous spectrum) in a ID double-barrier well potential 
in the absence of laser field. Details will be reported elsewhere [42]. With these 
solutions, we calculate the dipole matrix element, (a;)c-g, of the continuum-ground- 
state transition, which plays a central role in the resonant HHG emission [SI [201 ^ [TD] . 
The phase of this matrix element, shown in figure [5| jumps rapidly by 2.65 rad at 
the exact resonance position, Crcs = 26.27 eV. In both the 'detuned' and 'resonant' 
cases, H17 is still on the "blue" side of this jump, whereas H15 is on the "red" 
side. The recombination matrix element thus adds an additional phase difference 
of AipIq^ « 2.2 rad between H15 and H17, and thus shifts the GD at SB16 by 
AlpIq^ /2uji^ — 460 as. This is in very satisfactory agreement with the GD-shifts 
observed in the experiment and the numerical model. Note that the dispersion of the 
HHG medium around the resonance leads to a phase distortion with sharp excursions 
of opposite sign on either side of the resonance, qualitatively different from the step- 
like phase induced by (x)c-g. If this propagation-induced phase were important, the 
corresponding group delay shift at SB16 would have to be accompanied by a similar 
shift with opposite sign at SB18. Our experimental GDs are thus compatible with 
(i)c-g being the dominant source of phase distortions. 

The matrix element phase can also qualitatively explain the suppression of the 
RABBiT-oscillations in the 'resonant' case. As visible in figure[5) in the 'detuned' case, 
it varies very little within the width of H17, while in the 'resonant' case, the onset 
of the rapid phase jump causes a much stronger phase variation within the width 
of H17. This will contribute to the observed suppression of the RABBiT-oscillation 
contrast on the sidebands surrounding H17. The first order of the phase variation 
(linear component) corresponds to a delay of H17 with respect to the other harmonics. 
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For the 'resonant' case, the delay predicted by {x)c-g is 1.5 fs, which agrees well with 
the numerical TDSE solution. Unfortunately, this small value cannot be resolved in 
our experimental cross-correlations. With H17 exactly on resonance (Al = 802 nm), 
{x)c-g predicts a delay of 8 fs, which would certainly be measurable. 

Other effects could also contribute to the reduction of the contrast of the RABBIT- 
oscillation. For example, the resonant H17 wavefront could be spatially distorted, 
which would average out the relative phases of harmonic and probe beams in the mbes 
sensitive volume. Such a distortion, affecting only the resonant harmonic, could be 
caused by the dispersion around the resonance frequency if there were strong gradients 
of the Sn"'"-density in the hhg medium. But, as discussed above, the propagation- 
induced phase distortions seem negligible in our generation conditions due presumably 
to low plasma densities. 

4. Conclusion 

In conclusion, our study gives an affirmative answer to the practical question of 
whether resonance-enhanced HHG is indeed a source of intense ultrashort XUV pulses. 
The enhanced harmonic order has the same femtosecond duration as the non-resonant 
ones. On the attosecond time-scale however, we find significant distortions of the phase 
of the near-resonant harmonic. Our results suggest the detuning from the resonance 
as an effective handle controlling the resonant harmonic phase. 

From a more fundamental viewpoint, previous studies of the HHG phase 
properties focused mainly on the phase accumulated by the quasi-free electron in 
the continuum, or on the recombination step as a probe of molecular structure 
and dynamics |43l I44j . This article presents experimental evidence of the dramatic 
influence of the recombination step on the phase of resonant harmonics from an atomic 
target. 

We have shown that the recombination dipole matrix element alone can describe 
the origin of the phase distortions observed in our experiments. This confirms the 
four-step picture [20] illustrated in figure |4j where the resonance comes into play 
only as an intermediate state in recombination, but not as an initial state as in 
the model proposed in [22]. It is possible to factor out this recombination dipole 
matrix element from the expression for the radiating dipole in the description of 
HHG [U [To]. This is the basis for the so called 'self-probing' schemes to extract 
structural and dynamic information about the generating system from intensity, 
phase and polarization measurement of high harmonics [121 [Hj. Our results thus 
suggest the possibility of devising 'self-probing' schemes for atomic resonances based 
on the advanced characterization of resonance-enhanced high harmonics. In particular, 
the rapid phase variation responsible for the intriguing suppression of the RABITT 
oscillations may encode characteristic features of the involved AI state. To gain 
further insight, the phase within the spectral width of the resonant harmonic could be 
measured by XFROG [IS], Xuv-SPIDER [1^ or the FROG-CRAB [17] method, and the 
resonant harmonic could be spatially characterized [H]. 
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